Classical Heisenberg antiferromagnet on a garnet lattice: a Monte Carlo simulation 
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We have studied a classical antiferromagnet on a garnet lattice by means of Monte Carlo simula- 
tions in an attempt to examine the role of geometrical frustration in Gadolinium Gallium Garnet, 
Gd3GasOi2 (GGG). Low-temperature specific heat, magnetisation, susceptibility, the autocorrela- 
tion function A(t) and the neutron scattering function S(Q) have been calculated for several models 
including different types of magnetic interactions and with the presence of an external magnetic 
field applied along the principal symmetry axes. A model, which includes only nearest-neighbour 
exchange, Ji , neither orders down to the lowest temperature nor does it show any tendency towards 
forming a short-range coplanar spin structure. This model, however, does demonstrate a magnetic 
field induced ordering below T ~ O.OUi. In order to reproduce the experimentally observed prop- 
erties of GGG, the simulated model must include nearest neighbour exchange interactions and also 
dipolar forces. The presence of weak next-to-nearest exchange interactions is found to be insignif- 
icant. In zero field S(Q) exhibits diffuse magnetic scattering around positions in reciprocal space 
where antiferromagnetic Bragg peaks appear in an applied magnetic field. 



PACS numbers: 75.50.Ee, 75.40.Mg 
I. INTRODUCTION 



The introduction of frustration to magnetic systems 
leads to extra degeneracy for the ground state in addi- 
tion to the degeneracy resulting from the symmetry of 
magnetic Hamiltonian. The larger this additional de- 
generacy, the more likely frustration is to cause dramatic 
changes in the magnetic properties of the system, such as 
the absence of long range order even at the lowest tem- 
perature. Geometrical frustration has been one of the 
key issues in magnetism for atJeast twenty years. A re- 
cent wave of theoretical paperaMl as well as publications 
dealing with Monte Carlo simulationETD has emphasised 
the unusual magnetic properties of geometrically frus- 
trated systems. The question of whether the frustration 
leads to a disordered gapped state or to long-range Need 
type order in different types of geometry is still under 
debate. Current efforts seem to be concgoixated around 
two types of lattices: the Kagome lattice3'EH3 and the py- 
rochlore latticeffi. Recently has it been established that 
the pyrochlore lattice represents the only simple system 
for which the additional degree of freedom caused by the 
frustration is extensive - it is proportional to the number 
of spins involvecfi. 

The growth of theoretical interest in the pyrochlore 
lattice, a lattice of corner-sharing tetrahedra, is driven 
largely by experimental discoveries!! There are many 
chemically clean pyrochlores (some of which may be pro- 
duced as single crystalalTJ) with different types of mag- 
netic atoms and interactions, which allows one to pick 
the most suitable one for study and for comparison 
with a particular theoretical model. By studying the 
phenomenon in general a much better understanding of 
the magnetic properties of individual compounds can be 
achieved. The same reasoning applies to another geo- 
metrically frustrated system - an antiferromagnet on a 



Kagome lattice, where SrCrapGai2-9 P Oig0, jarosites@ 
and some other compoundsli3 provide quite a variety of 
model systems. 

Gadolinium gallium garnet, Gd3Ga50i2, is a unique 
example of an antiferromagnet on the garnet lattice. 
There are no other compounds matching its magnetic 
properties. In GGG (space group Ia3d) the magnetic Gd 
ions are located on two interpenetrating, corner-sharing 
triangular sublattices, where the triangles of spins do not 
lie in the same plane - the angle between two nearest tri- 
angles is equal to the angle between the diagonals of a 
cube, 70.5° (see Fig. Q). In this compound the trian- 
gular arrangement of the nearest spins is combined with 
complete exchange isotropy (the single-ion anisotropy is 
negligibly smalta) and with a relatively strong dipole- 
dipole energy. Although the magnetic properties of vari- 
ous garnets have been thoroughly studied during the past 
half century, the analogy between any of them and GGG 
is not straightforward. All other magnetic garnets order 
at some low temperature, while GGG does not. No long 
range magnetic order has been detected in GGG down to 
25 mKEj, while other gallium garnets based on Dy, Nd, 
Sm and Er, rather than Gd, have been found, to be mag- 
netically ordered at temperatures below 1 KLLj. The near- 
est analogy to GGG would probably be found among the 
Mn-based garnets, where the single-ion anisotropy is also 
very small. However two similar magnetically isotropic 
garnets, Mn3Al2Ge30i2 and MnaA^SisOi^EZI, also order. 
Most likely this is due to the presence of relatively strong 
next-to-nearest exchange interactions. If and when the 
degeneracy of the ground state is removed and the sys- 
tem undergoes a phase transition to a long-range ordered 
state, almost all complications disappear. The magnetic 
ground state and the main interactions are known from 
experiment and theoretical calculations are straightfor- 
ward. Numerical estimates exist to at least the accuracy 
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that experiments currently attain. However a theoretical 
model describing adequately the magnetic properties of 
GGG still has to be developed. 

This paper presents the results of classical Monte Carlo 
simulations of the magnetic properties of the Heisenberg 
antifcrromagnet on a garnet lattice. While some of the 
initial results related to the GGG have-been briefly re- 
ported in our neutron scattering papersEjO, where they 
have been used to explain the obtained experimental data 
and also to predict possible experiments, here we take a 
more general approach to the problem. We address is- 
sues which are a not necessarily directly related to GGG, 
but are interesting from a theoretical point of view, e.g. 
we discuss properties of a model which includes nearest- 
neighbour exchange interactions only. Where possible we 
compare with the results of simulations for the pyrochlore 
and Kagome lattices and show, that an antiferromagnet 
on a garnet lattice is yet another highly frustrated mag- 
netic system exhibiting a number of unusual and intrigu- 
ing properties. 



II. SIMULATION MODELS 



We consider the Hamiltonian 



n 



]T JsSiSj +d J2 



<hJ> 



<hJ> 



SiSj _ 3 (S i r ij )(S j r ij ) 



(1) 



where the spins Si are classical, three-component vectors 
on the Gd 3+ sites of a garnet lattice, S=7/2 as in GGG. 
The first term is the exchange interaction, the second 
term is the dipole-dipole interaction between the mag- 
netic moments. 

The original idea to simulate the magnetic proper- 
ties p£ GGG using MC methods belongs to Kinney and 
WoliE], who calculated the temperature dependence of 
the specific heat and by comparing the results with the 
experimental data have obtained the amplitudes of the 
nearest and next to nearest neighbour exchange interac- 
tions Ji, J r 2 and J3. More recently Schiffer et. al cal- 
culated the magnetic phase boundary and have inves- 
tigated the magnetic structure of GGG in an applied 
fieldEll. We use the same value of the nearest exchange 
constant as Kinney and Woha, J\ = 0.107 K, because it 
produces good estimates for the temperature dependence 
of the susceptibility and also for the saturation field of 
the magnetisatior£3. However, as will be shown later, the 
values of J 2 and J3 quoted in Ref£3 are not essential: as 
long as they are small in comparison with J\ , they do not 
change significantly the predicted magnetic properties of 
GGG and therefore can not be reliably determined from 
the MC simulations. 

The strength of the dipole-dipole interaction, D, is de- 
fined by the distances between the ith and jth spins. In 
GGG the Gd 3+ sites are separated by = 3.781 A, 



where a = 12.349 A is the lattice constant at low temper- 
ature. Therefore we use Ddd = 0.0457 K for the strength 
of the nearest-neighbour dipolar interaction. It is very 
important and at the same time very difficult to simulate 
reliably such a long range interaction as the dipole-dipole 
one. Eor some simpler lattices, for example, a 2D-square 
latticeEa, or forJaighly anisotropical systems, such as spin 
ice pyrochloreso, the Ewald summation technique can be 
used to treat the long-range nature of the dipole-dipole 
interaction. In case of Heisenberg spins located on the 
complicated lattice of GGG, however, there is no option 
but to introduce a cut-off range, i?o, and to neglect the 
dipole-dipole interaction for all distances larger than Rq. 
Previous simulations have restricted the dipole-dipole in- 
teraction to a third neighbouiErtil, while in our model Rq 
has been extended to include the fourth neighbour. We 
have also made several test runs to compare the sim- 
ulation results for this model with both shorter (to a 
third neighbour) and also longer (ten neighbours) cut-off 
ranges and have found no significant difference, which 
suggests that this model describes the dipolar force rea- 
sonably well. The dipolar interaction between two mag- 
netic moments decays as the number of neighbours 
in a shell SR is proportional to R 2 , therefore the dipole- 
dipole energy should decay only relatively slowly, as j=r. 
In reality, however, the extension of the cut-off range 
from a third neighbour to a fourth one does not change 
significantly either the total dipolar energy, nor the over- 
all system energy. A possible answer to this puzzle might 
be related to the fact that all magnetic interactions in 
GGG, including the dipolar one, are frustrated: the con- 
tribution of the individual magnetic moments to the to- 
tal system energy is mutually cancelled or nearly can- 
celled, therefore for each magnetic moment only the lo- 
cal surroundings influence the choice of magnetic orienta- 
tion. Similar observations have been made during recent 
Monte Carlo simulations on pyrochlore lattice which in- 
cluded long-ranged dipole-dipole as well as short-ranged 
exchange interactions!! 

MC simulations have been performed for lattice sizes 
Lx Lx L, with L — 3 to 9 unit cells, containing 648 to 
17496 spins. Significantly larger lattice sizes, than pre- 
viously used, have ensured that the magnetic correlation 
length in the disordered phase does not exceed the system 
size. Simulations with larger lattice sizes have improved 
the resolution of the calculated scattering function, S(Q), 
in an applied magnetic field allowing us to resolve clearly 
individual magnetic Bragg peaks. A standard Metropolis 
algorithm with periodic or open boundary condition has 
been employed; up to several millions Monte Carlo steps 
per spin (MCS) were performed at the lowest tempera- 
tures. Where possible an attenuation factor SS has been 
introduced in such a way that roughly 50% of the at- 
tempted spin moves were acceptedEj, which has resulted 
in a dramatic increase of the spin relaxation rate. For 
the simulations in a magnetic field this procedure has 
been abandoned to permit the system to make abrupt 



2 



structural changes. 

The magnetic field is assumed to be applied along the 
(100) direction unless otherwise stated. 



III. RESULTS AND DISCUSSION 



A. Zero external field properties 

We begin by addressing the issue of the phase transi- 
tion at low temperature in zero magnetic field. In GGG 
no sign of long i*ange magnetic order has been found 
down to 25 mKEj, moreover, frustration induced spin 
freezing has been suggested at temperatures below 125 - 
135 mK on the basis of single crystal magnetisation mea- 
surements: the susceptibility is frequency dependent, and 
the static magnetisation is different for field cooling and 
zero field coolingEj. However, neutron scattering experi- 
ments show that at the lowest temperatures the magnetic 
system is not frozen completelyta. It rather behaves as a 
mixture of a liquid and solid states. 

The first thing to notice is that the simulation model, 
which includes only nearest-neighbour exchange, J\ , does 
not show any sign of a phase transition down to at least 
T=l mK (which is less than 0.1% of the exchange energy 
JS 2 ). Several measured quantities show thai-jthe system 
remains in a spin-liquid (or, following VillainE3, a cooper- 
ative paramagnet) phase: averaging over sufficiently long 
time gives zero magnetic moment on each site, the scat- 
tering function S(Q) does not show any sharp peaks, the 
magnetic correlation length Q(r) =< S(0)-S(r) > does 
not exceed the system size (see Fig. ^). In fact, close in- 
spection of Fig. p| reveals that correlations are very small 
beyond the first unit cell. In addition there is no obvious 
maximum or cusp in the heat capacity temperature de- 
pendence (see Fig. |^). To test the suspicion that at low 
temperature Monte Carlo simulations are not effective 
enough in allowing the system to reach equilibrium, we 
have checked whether the simulation results depend upon 
the starting conditions. No difference in results have been 
noticed when starting calculations from an initially ran- 
dom or a 120°-degree planar triangular state. 

The low-temperature specific heat itself is an impor- 
tant thermodynamic quantity, whose value is sensitive to 
the presence of zero modesa and quartic modest! In the 
pyrochlore lattice each quadratic mode contributes fcs/2 
to the heat capacity, each quartic mode fce/4 and zero 
modes do not contribute at all, thus, reducing the zero- 
temperature specific heat to 3fce,/4El, while in Kagome 
lattice it is reduced to llfcs/1219. Our initial calcula- 
tion on a relatively small system with periodic bound- 
ary conditions showed that C(T =0) was indistinguish- 
able from unity within the accuracy of the simulations. 
However, prompted by the comparison with the Kagome 
and pyrochlore lattice results, we have performed much 
longer Monte Carlo runs on much bigger systems with 
open boundary conditions (we use open boundary con- 



ditions in order to avoid the imposition of periodicity on 
a potentially incommensurate magnetic system). As can 
be seen from Fig. |, C(T = 0) « 0.94(2) with the accu- 
racy of the calculation sufficiently high to claim that it 
is actually below unity. There is no significant difference 
in C(T =0) calculated for systems of 5x5x5 and 9x9x9 
containing 3000 and 17496 spins respectively. 

The introduction of the dipolar interactions slows 
down the spin-relaxation process. Fig. displays the 
time dependence (time is measured in MCS) of the auto- 
correlation function A(t) — -h < Sj(0)Sj(t) > for the 
two models: with (bottom) and without (top) dipolar 
forces. The model which includes dipole-dipole interac- 
tions does not show noticeable relaxation by T = 50 mK, 
while the model with only nearest exchange interactions 
is still relaxing even at an order of magnitude lower tem- 
perature. The difference between the autocorrelation 
function for these two models is evident at all temper- 
atures below 0.5 K, which approximately coincides with 
the nearest neighbour dipolar energy, D^d X S 2 . 

Dealing with very slow relaxing spin systems and a po- 
tential spin-glass transition it is essential to ensure that 
the simulation time is longer than the equilibration time. 
In practice the first to MCS are used only for equilibra- 
tion and then calculations and averaging are carried out 
during the next to steps. An estimation of an appropri- 
ate value of to could be obtaincdJbllowing the procedure 
introduced by Bhatt and Youngc3, where the spin-glass 
susceptibility, XSG ; has been calculated in two different 
ways. In the first method we calculate an overlap be- 
tween two uncorrelated sets of spins which approaches 
XSG from below, if to is shorter than the equilibration 
time. In the second approach the four-spin-correlation 
function is calculated, which approaches XSG from above, 
if to is small. The to is considered to be long enough 
and the results arc accepted only if the two estimates of 
Xsg agreed. In a GGG model which includes both ex- 
change and dipolar interactions to becomes enormously 
long at low temperatures. In fact even during the runs 
with t = 10 6 MCS the results showed no agreement 
between the two approaches for all temperatures below 
T = 100 mK. Therefore the results of calculations in zero 
field for a model which includes dipole-dipole interactions 
could not be considered as reliable below this tempera- 
ture. The problem of long equilibration times is removed 
by the application of an external magnetic field. 

The results of the simulations with the model, which 
takes into account only the nearest neighbour exchange 
interaction, J\ , fits well the experimental neutron scattet- 
ing function S P (Q) at all temperatures above 140 mKH. 
The neutron scattering function S P (Q) for a powder sam- 
ple is calculated as: 



S P (Q) = f{Q) 2 Y.< s * s i 

i,3 



> 



sm(Qr) 
Qr ' 



(2) 



where f(Q) is the magnetic form factor. S P (Q) has sev- 
eral broad diffuse scattering peaks (see Fig. 5 in Ref.lla) , 
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whose intensity increases as the temperature decreases in 
agreement with the experiment. The introduction of the 
dipole-dipole interaction at these temperatures does not 
change S P (Q) significantly. 

A somewhat unexpected results have been obtained 
earlieicj for a single crystal neutron scattering function, 
calculated as: 

N 

S x t(Q) = (f(Q)J2^e iQr -) 2 , (3) 

n 

where q„ is the magnetic interaction vector. Even at 
temperature well above T = 140 mK, where there is 
no problem from very long equilibration times, S x t{Q) 
demonstrates a tendency to form incommensurate peaks 
around intege*-i positions in the reciprocal space (see 
Fig. 4 in Ref.Ej). The intensity of these incommensu- 
rate peaks is much lower than the expected intensity of 
the true long-range order Bragg peaks, and their width 
is determined by the system size. The exact position of 
these peaks in reciprocal space is not fixed, it may change 
from one "snap shot" of S x t(Q) to another. Only after 
averaging significantly large amount of the "snap shots" 
(from several dozens to several hundreds) a clear picture 
of the short-range incommensurate magnetic order was 
obtained. However, this is most likely to be an artifi- 
cial effect caused by the periodic boundary conditions: 
when they are removed, the effect of splitting seems to 
disappear. Top panel of Figure ^| shows simulated single 
crystal neutron scattering intensity of GGG in the (hkO) 
plane at T = 0.2 K. 

Another interesting aspect of this study is to investi- 
gate how the ratio of exchange to dipolar interactions in- 
fluences properties of the Heisenberg antiferromagnet on 
a garnet lattice. In GGG J\ is about twice the strength 
of Ddd OTd there is no magnetic order, while in Mn-based 
garnetsEZl the ratio J\ / Ddd is slightly higher and they do 
order. For instance, in Mn3Al2Ge30i2, which undergoes 
an antiferromagnetic phase transition to a 120°-structure 
at Tn = 6.65 K, the J\ =0.57 KEj is more than ten times 
stronger than Ddd in GGG. Our simulation shows that 
this fact alone could not lead to the appearance of the 
long-range magnetic order. In a model, where J\ has 
been increased up to a hundred times keeping the Ddd 
value fixed, the ground state remained disordered. How- 
ever, the introduction of the nepdj-to-nearest exchange 
interaction with a value cited ir£3, J 2 — 0.12 K, does 
make a difference: the system immediately undergoes 
a phase transition to a LRO state, which reveals itself 
clearly both as a cusp in a heat capacity temperature de- 
pendence and as peaks in the scattering function, S(Q). 

B. Magnetic properties in an applied field 

As has been mentioned above, in an applied magnetic 
field the problem of long equilibration times is much 
less severe, which gives us an excellent opportunity to 



investigate the magnetic phase diagram of GGG in de- 
tail. A phase transition to a LRO state in magnetic field 
was detected by calculating the specific heat tempera- 
ture dependence in constant field or by calculating its 
field dependence at constant temperature. Figures 4 and 
5 in Ref.EJ give examples of such calculations. The po- 
sition of the specific heat maximum is not sensitive to 
the introduction of the relatively weak next-to-nearest 
exchange interactions (such as were quoted in Ref.EJ, 
J 2 = -0.003 K and J 3 = 0.010 K), neither does it show 
any visible size-dependence. In the field dependence of 
the specific heat, only one anomaly corresponding to the 
upper transition field is well-pronounced, while there is 
no obvious anomaly corresponding to the lower transition 
field, which agrees with previous MC simulationsEU. 

In order to reproduce accurately the experimentally 
observed phase diagram of GGG, the simulation model 
must include nearest neighbour exchange interactions 
and also dipolar forces. However, even in the model in- 
cluding only nearest neighbour exchange more accurate 
calculations revealed signs of the phase transitions in a 
magnetic field. Fig. ^ presents the magnetisation curves 
and also their derivatives at T = 1 mK for such mod- 
els. Before reaching a saturation point at if « 1.7 T, the 
raw magnetisation shows a relatively small change of the 
slope around if » 0.6 T, which is not really a striking 
feature and therefore has passed unnoticed in our ear- 
lier calculations. In the susceptibility curves, however, a 
clear minimum is present at if ~ 0.6 T. We believe that 
this minimum in susceptibility corresponds to the ap- 
pearance of a collinear long-range ordered state induced 
by an applied magnetic field. In complete agreement with 
the theorjO, which analyses an order by disorder mech- 
anism in various highly frustrated antiferromagnets, an 
ordering happens only around a special value of the mag- 
netic field - one third of the saturation field in case of a 
garnet lattice. 

In GGG, that is in a model which includes the dipolar 
forces, an ordered magnetic structure induced by an ap- 
plied field is characterised by the appearance of a nonzero 
average value of the perpendicular component of local 
magnetisation. The field dependence of parallel and per- 
pendicular components of local magnetisation at con- 
stant temperature is shown on Fig. f^, while Fig. |8| shows 
its temperature dependence in constant field. The two 
sets of curves on each of these figures reflect the fact that 
in applied magnetic field the 24 Gd sites are split un- 
equally into two different syrinmetry sites - "A" and "B" 
sites in the notation of Rcf.Eil. When the field is applied 
along the (001) direction, the 8 "A" sites (represented 
by solid symbols on Fig. [?] and ^) are of higher symme- 
try than the 16 "B" sites (represented by open symbols). 
Clearly, < M > xy on the "A" sites serves as on order 
parameter for the transition from a paramagnetic state 
into an antiferromagnetically ordered state. Spins on the 
"B" sites, however, retain a nonzero value of the perpen- 
dicular component of magnetisation even in the param- 
agnetic state. This effect is caused by the dipole-dipole 
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interaction. 

In order to avoid problems with possibly many 
metastable states the calculations were always started 
at high temperatures and fields and then the system an- 
nealed as it came into equilibrium at the desired field 
and temperature for measurement. However, even tak- 
ing these precautions the problem of long equilibration 
times at low-temperature low-field region was unavoid- 
able. Therefore an abrupt jump of magnetisation around 
H = 0.25 T clearly visible on Fig. ^ is most likely to be 
an artificial result. 

Even without detailed knowledge of the magnetic 
structure in a field we can check how stable it is to the in- 
troduction of second and third next to nearest exchange 
interactions, J2 and J3. This has been done by calculat- 
ing the neutron scattering function S P (Q) for a 3x3x3 
system according to formula (||) . The results suggest that 
the magnetic order is rather stable in all four quadrants 
in the J2 — J3 plane. The structure does not change 
when J 2 = -0.003 K and J 3 = 0.010 K axe introduced 
corresponding to the values quoted in Refo. 

In an applied magnetic field, where LRO is developed, 
the formula (^) is no longer valid. Although it unambigu- 
ously shows the appearance of magnetic Bragg peaks, 
their intensity is not calculated correctly. However the 
overall field dependence of tlic.intcnsity mimics extremely 
well the experimental datata There are two different 
groups of magnetic peaks, ferromagnetic and antiferro- 
magnetic. The intensity of the former group is growing 
in lower fields and saturating in a higher field. The inten- 
sity of the latter group also grows in lower fields reaching 
a maximum at around H = 1 T and then decreases in 
higher fields and disappears above H = 2 T. Exactly 
the same behaviour has been seen by simulating single- 
crystal scattering intensity according to formula (^). Bot- 
tom panel of Fig. || shows the results of such calculations 
for H — 1.06 T and T = 200 mK. In an applied mag- 
netic field a set of strong and sharp Bragg peaks replaces 
diffuse magnetic scattering observed in zero field. One 
interesting aspect of the calculations must be empha- 
sised here: the relative intensity of the symmetry related 
antifcrromagnetic peaks, such as, for example, (210) and 
(120), is not constant in time. The intensity of each peak 
may change arbitrarily at any time from almost zero up 
to maximum value, while the sum of two peaks intensities 
remains constant. 

The only discrepancy between the Monte Carlo results 
and the neutron scattering data in magnetic field is the 
presence in the later of an incommensurate peak located 
between two antiferromagnetic peaks, (200) and (210). 
While this incommensurate peak is clearly visible in the 
neutron scattering dataEfl, neither S P (Q) nor S x t{Q) (in- 
cluding a model where the magnetic field is applied along 
the (110) and (111) directions) demonstrates peaks at an 
incommensurate position. The reason for the discrepancy 
remains unknown at the moment. 



IV. CONCLUSIONS 

To summarise, we have presented the results of clas- 
sical Monte Carlo simulations for the low temperature 
behaviour of the frustrated antiferromagnet on a garnet 
lattice. We have studied several different models, pay- 
ing particular attention to two of them. The first model, 
which includes only nearest-neighbour exchange interac- 
tions, does not order down to lowest temperature, neither 
does it show any signs of spin-freezing. Calculations of 
the zero-temperature specific heat for such model suggest 
the presence of soft modes. The indications of the phase 
transition into an ordered (presumably collinear) state 
have been found at low temperature in applied magnetic 
field around a third of the saturation field. The exper- 
imentally measured properties of GGG including the H 
vs. T magnetic phase diagram are consistent with our 
findings for the simulation model which includes nearest 
neighbour exchange interactions and also dipolar forces. 
A perpendicular component of local magnetisation serves 
as an order parameter for the phase transition in an ap- 
plied field. 

In conclusion we discuss several questions, which have 
been considered in this article, but which most certainly 
require further theoretical investigations. 

Firstly, both the low-temperature specific heat, C(T = 
0), and the single crystal scattering function, S x t(Q), are 
unusually sensitive to the boundary conditions. The in- 
fluence of particular boundary conditions on the appear- 
ance of soft modes and incommensurate peaks in S x t(Q) 
needs to be examined further. 

Secondly, in the ordered state only the total intensity 
for the pairs of a symmetry related antiferromagnetic 
Bragg peaks, remains constant, while the intensity of an 
individual peaks may change. Could this behaviour be 
related to the energy-free motion of long chains of mag- 
netic moments, similar to what happens in the Kagome 
lattice, or is it a sign of the domain walls movement? 

We are grateful to J. T. Chalker, M. J. P. Gingras, 
M. J. Harris, C. L. Henley, P. C. Holdsworth, R. Moess- 
ner and M. E. Zhitomirsky for discussions and other valu- 
able contributions and also to J. N. Reimers for his help 
with the computer calculations on the initial stages of 
the project. 
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FIG. 1. Positions of the magnetic Gd ions in a garnet struc- 
ture. There are 24 magnetic ions per unit cell, they are di- 
vided into two interpenetrating sublattices. 
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FIG. 2. Correlation function for a system of 5x5x5 unit cell 
sizes (3000 spins), which includes only nearest neighbour ex- 
change interaction, Ji. Top picture shows correlation between 
spins belonging to the same sublattice, bottom - correlation 
between sublattices. 
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FIG. 4. Time dependence (in Monte Carlo steps per spin) 
of the autocorrelation function < S(0)S(t) > for various 
temperature from T = 1 K down to 2 mK in a model which 
includes: a) only the nearest neighbour exchange interaction, 
Ji, b) the nearest neighbour exchange interaction and the 
dipole-dipole interactions up to fourth neighbour. System 
with a lattice size 9x9x9 unit cells has been used for this 
calculations. 
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FIG. 5. Simulated single crystal neutron scattering inten- 
sity of in the (hkO) plane at T — 0.2 K in a zero field (top 
panel) and in a field of H — 1.06 T applied along (001) direc- 
tion (bottom panel). The data have been obtained according 
to formula ^ for a model size of 9x9x9 unit cells with open 
boundary conditions. 




FIG. 6. Field dependence of the magnetisation (top) and 
susceptibility (bottom) at T = 2 mK, H \\ (001). A 5x5x5 
model has been used to generate these data. Open and solid 
symbols represent the data for the model which included near- 
est neighbour exchange interaction with and without dipolar 
forces respectively. 
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FIG. 7. Field dependence of parallel (bottom) and perpen- 
dicular (top) component of local magnetisation at T—0.1 K, 
H || (001). For notation see main text. 
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FIG. 8. Temperature dependence of parallel (bottom) and 
perpendicular (top) component of local magnetisation in a 
field # = 1.06 T, H \\ (001). For notation see main text. 
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